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Abstract 

We present a new class of flat-band Hubbard models which have saturated ferromag- 
netic ground states at two distinct electron numbers for different values of parameters. 
The models are extensions of Tasaki's flat-band models. 

1 Introduction 

It is widely believed that the spin-independent Coulomb interaction and the Pauli exclusion 
principle can generate ferromagnetism in itinerant electron systems. One of the motivations 
to study the Hubbard model has been to establish and understand the generation of ferro- 
magnetism in simplified situations taking account of these effects p. Mielke and Tasaki 
jHj, independently, presented Hubbard models which exhibit saturated ferromagnetism for 
certain electron numbers when the Coulomb interaction U is positive. These models have 
a common feature that the single electron spectra contain dispersionless bands, and are 
called fiat-band models. In j4j and 0, Tasaki also discovered Hubbard models exhibiting 
ferromagnetism which models are nonsingular in the sense that both the density of states 
and the Coulomb interaction are finite. Recently Tanaka and Ueda succeeded in proving 
the existence of saturated ferromagnetism in a Hubbard model obtained by adding extra 
hopping terms to Mielke's fiat-band model on the kagome lattice 

Although the fiat band Hubbard models are singular and not physically realistic, their 
study can be a basis of more realistic results about ferromagnetism. It is therefore important 
to find out which flat-band models exhibit ferromagnetism. Although an abstract criterion 
was presented by Mielke j7j, we still do not know precise class of models which satisfy the 
criterion. 

In this paper, we follow Tasaki's construction of his flat band models, and construct a new 
class of Hubbard models in arbitrary dimensions with finite U and finite-range hopping. We 
prove that the models exhibit ferromagnetism in their ground states at two distinct electron 
numbers. The difference between Tasaki's original model and ours can easily be seen from 
Fig^ where the simplest one-dimensional versions of the models are illustrated. Tasaki's 
model has one "internal site" (gray dot in the figure) in each unit cell, while ours has two. 
This difference in lattice structure makes our model to have different "exchange mechanism" 
where a single-electron state localized at each pair of internal sites play an important role. 
See section HI for more details. 
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Figure 1: The lattice structure and the hopping amphtude in the one dimensional flat band 
models of (a) Tasaki's and (b) ours. The black dots are the external sites (in S) and the 
gray dots are the internal sites (in X). Tasaki's model has one internal cite in each cell while 
our has two. 

We have thus found that a extension of Tasaki's construction lead to a new class of 
models exhibiting ferromagnetism. We hope this study will shed light on general structure 
of fiat-band ferromagnetism. 



2 The model and main results 
2.1 Construction of the lattice 

In the original fiat-band models by Tasaki the basic cell in the lattice A consists of a 
single internal site and some external sites. In our new models, the basic cell consists of two 
internal sites and some external sites. More precisely we let the basic cell be 

C = {u,V,Xi,X2,- ■ ■ ,Xn}. (2.1) 

We call u and v the internal sites of C, and Xi,X2, - ■ ■ ,Xn the external sites. 

To form the lattice A, we assemble M identical copies Ci, C2, • ■ ■ , Cm of the basic cell C, 
and identify external sites from m distinct cells and regard them as a single site. In other 
words, an external site in A is shared by m distinct cells. We denote by |A| the number of 
sites in A. See FigI21for an example of a cell and a resulting lattice. 

The lattice is naturally decomposed as 

A=IUS (2.2) 

where I and S are the sets of internal sites and external sites, respectively. We also denote 
by JT" the assembly {1, 2, ■ ■ ■ , M} of the indices of cells. From the above construction, we see 
that the numbers of sites in these sublattices are \I\ = 2M, \S\ = nM/m. By using |A|, we 
can write |X| = 2m|A|/(2m + n), \S\ = n\A\/{2m + n) and = m|A|/(2m + n). 
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In what follows we always regard Cj as a subset of A. We denote the two internal sites 
in Cj as Uj and Vj. For an external site x G we denote by the collection of indices j 
such that X & Cj. We also define C A to be the union of m cells which contain the site x. 




Figure 2: An example of a cell and a lattice. From the quadrangular cell with four sites 
(a), one can form (b) a decorated square lattice by identifying four external sites. This 
defines flat-band Hubbard model in two dimensions. For the precise definition see (I2.19|) 
and (|2.2Up (where one should set ri = 2, m = 4). When = |A|/5, the model exhibits 
saturated ferromagnetism for any t > 0, s' > and U > 0. The model also exhibits saturated 
ferro magnetism for any t > 0, U > when A^^e = 3|A|/5 and s' = 0. See Theorem 12.11 and 
Theorem 12.21 

2.2 Fermion operators 

We consider an electron system on the lattice A. For each site r G A and a =|, |, we define 
the creation and the annihilation operators cj^ and Cr^a for an electron at site r with spin a. 
These operators satisfy the canonical anticommutation relations 

{'-'r,(T' ~ ^r,s^(7,Ti (2-3) 

and 

{4,<x,<r} = K.,c.,.} = 0, (2.4) 

for any r, s G A and a, t =|, |, where {A, B} = AB + BA. We denote by $vac a normalized 
vector state which satisfies Cr,cr*^' vac = for any r G A and a =t, i- Then for arbitrary subsets 
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A|,A| C A, we define a state 



n n *. 



(2.5) 



, rGA-r 



, rGAi 



in which sites in A| are occupied by up-spin electrons and sites in A^ by down-spin electrons. 
Next we define total spin operator Stot = ('S'tot; "^tot' '^tot) by 



^tot — r, 2^ ^r,a\P )(j,tW,t 



(2.6) 



for a = 1,2, and 3. Here p*-"-* are the Pauli matrices defined by 

^"^(°o)' 

We finally define special fermion operators as in [Hj. Let > be a constant. For x & S, 

(2.8) 



let 



where the sum is over m sites adjacent to x. For j G JT", let 



(2.9) 



where the sum is over the n external sites in the cell Cj, and 

dj^a C^i^ (J Ci]j^(j. (2.10) 

From the anticommutation relations for the basic c operators, one can easily verify that 

h,r} = {4,., d^A = 4,.} = (2.11) 

for any x E S, j, k E J ^ and a, r =t, |. 

The anticommutation relations for the a operators are 



mv 



1 H , if X = cr = r; 



lo. 



if X 7^ y, a = r; 
otherwise. 



(2.12) 



For x,y e S, we defined 

= lAxHA^nxi, (2.13) 

which is the number of the internal sites directly connected both to x and y. For the b 
operators, we similarly have 

{2 + nu'^, \l j — k, a — r; 
iij^k.a^T- (2.14) 
0, otherwise. 

For j, k e J, we defined 

ej^k^lCjOCkDS], (2.15) 

which is the number of external sites which are included in both Cj and C^. For the d 
operators, we have 

{4..4.} = {o: ottewtr^'^ (2.16) 

Since the states aj.^^$vac) &],o-^vac said (ij_^$vac are linearly independent and the number 
of these states are 2|A|, An arbitrary many-electron state of the system can be represented 
as a linear combination of the basis states 

^o(-Et, El, J|, Ji, Jp J[) = 

u<] (u<] Iw'i] Im.] Im] (n40*™c(2.i7) 

with arbitrary subsets i^i, C J|, Jj^, J^, J| C JT". Here J|| + | Jj^| + | J:[|-|-| | — 

Ng is the total electron number. 

2.3 An extension of Tasaki's flat-band model 

We study a Hubbard model with the Hamiltonian 

M M 

^ = ^ E E ^Ih'^ + ^' E E 4,^,^^ +uj2 <T^ni (2-18) 

where t > 0, s' and U > are real , and rir^a = cJo-'^r.o- is the number operator. We can 
rewrite the same Hamiltonian in the more standard form as 

^ = E tr,sCl,aCs,a + ^ E ^nT^^-i (2-19) 
r,seA reA 



where tr,s are the hopping amphtudes given by 



f t 

t 
t 



Ui,Ui 



Ui,Vi 



^X,Ui 



t + s', 
t - s', 



tx,y ~ 



tui,x tx,Vi 



t 



Vi,X 



0, 



0, 



X eCi 
X ^ Ci 



xe£ 
i e J 
i e J 

X e £ ,i e J] 

x,y eS,x^y; 



(2.20) 



From the anticommutation relations ()2.11|) and ()2.16p . one easily find that the single- 
electron Schrodinger equation corresponding to ()2.18j) has |£^|-fold degenerate eigenstates 
aj. ^.$vac with energy and M-fold degenerate eigenstates (ij^^$vac with energy 2s'. It is 
also easy to see that the remaining eigenvalues are positive. We have thus defined a new 
class of flat-band Hubbard models. The models exhibit ferromagnetism as the following two 
theorems state. 

Theorem 2.1 Consider the above Hubbard model with Nc = \S\ = n\A\/{2m + n) and 
s' > 0. For any U > 0, the ground states have total spin Stot = 'S'max(= Nc/2), and are 
non- degenerate apart from the trivial (2S'max + l)-fold degeneracy. 



Theorem 2.2 Consider the above Hubbard model with = \S\ + \J'\ = 
and s' = 0. For any U > 0, the ground states have total spin S^Q^ — 
non- degenerate apart from the trivial (2S'max + ^)-fold degeneracy. 



(n+m)|A|/(2m + n) 
x(= Nc/2), and are 



It is remarkable that the new models show saturated ferromagnetism at two distinct 
electron numbers for different values of the parameters. This is a unique property of our 
models. 



3 Proof 

We define the states $2t as 



\xee / 



2T 



n 



i<] n4t)*- 

\x&£ / \j&J 



We decompose the Hamiltonian as H = iJhop + -f^int where 

M M 

i=i 'T=r,i j=i a=r,i 



(3.1) 



(3.2) 
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rGA 

Note that both ifhop and ifint are positive semidefinite. 



(3.3) 



3.1 Proof of Theorem [O 

We consider the case with N^, = \S\ and s' > 0, and prove Theorem 12.11 Since the proof is 
essentially the same as that found in jHEj, we shall be brief. Since H >0 and if\E'i^ = 0, 
we see that an arbitrary ground state $gs satisfies 

i^hop<f GS = 0, i^int^GS = 0. (3.4) 

From the second relation in ()3.4|) . we further find that $gs must satisfy 

Cr,^Cr,l^GS = (3.5) 

for any rGA. 

By using ()2.17j] and the first condition in ()3.4j] . $gs can be represented as a linear 
combination of the basis states 

^,iE^,E^,J^,J^) = I n <T I n <i (3.6) 
By using the anticommutation relations {cx,a, (^LA = '^o-.r'^x.j/, we see that 



for any x E S. By using ()3.5|) for r G £^ we find that only the basis states satisfying 
E-^nEi =0, contribute to $gs- 

In this way, $gs can be written as 

$GS = E^M(n«U))'^- (3.8) 
a- \xee / 

where the sum is over all the spin configuration cr = {(Tx)xe£ on S and g[cr] is a coefficient. 

By using ()3.8|) and the anticommutation relations {cu^^a, ct^.r) = —{^/'^)^ct,tX[^ E S H Cj], 
for any x E £, where x[true] = 1, x [false] = 0, we get 

%,T%,i*GS = ^ sgn[a,/3](^[o-] -^[o-^^^]) j JJ «1,<t(x) I "^va^S.g) 

s.t.a>l3 s.t.<T(a)=t,<T(/3)=i 



where we have introduced an arbitrary ordering in 8 to avoid double counting and the factor 
sgn[a,/3] comes from the exchange of fermion operators. The spin configuration cTq,^/? is 
obtained from cr = {aa)aee by switching cTq, and ap. Since the basis states in ()3.9|1 are all 
linearly independent, we find from the property ()3.5|1 that 

g[(T]=g[(T^^p] (3.10) 
for the sites which belong to £r\Cj. Since the entire lattice is connected, ()3.10p ensures 

(3) 

that the lowest energy is unique in each sector with a fixed Sl^l- Therefore $11- is the unique 
ground state apart from the degeneracy for rotational invariance. This completes the proof 
of Theorem 12.11 



3.2 Proof of Theorem 12721 

We treat the case with = \S\ + \ J'\ and s' = 0, and prove Theorem 12.21 By using $2T 
instead of we find that the conditions ()3.4|) and ()3.5p are still valid. By using ()2.17|) and 
the first condition in ()3.4|1 . we find that an arbitrary ground state $gs can be represented 
as a linear combination of the basis states 

^,{E,, E,, J,, = ( n 4,T I n <i ( n 4t ( n 4,1 -^vac. (3.11) 

As before, (|3.5p and (j3.7p imply that only the basis states with = contribute to 

'^'gs- We will also prove that the basis states should satisfy J| fl = in order to contribute 
to $Gs- In other words, d states cannot be doubly occupied in a ground state. We prove 
Theorem 12.21 assuming this claim. 

From 1)3.111) and the above mentioned constraints, we have 

$GS = fn4,.(.)) f n ^vac (3.12) 

<T \xee / \jej / 

where the sum is over all the spin configurations cr = {crr)r£euj and g[<T] is a coefficient. Here, 
for notational simplicity, we identified the index set J' with the set of sites {uj}j(zj. 

By using ()3.12p and the anticommutation relations {c„^,cr, al ^} = — (z//2)5o-_^x[a; E £ nC 



0\i 



{cuj,a, 4,r} = ^^,rX[uk ^ Cj] for auy X e£, j,k e J, we get 

c«„TC«„i'^GS = Y sgn(a,/3)(^[o-] -^[(T^^^j) 



a,f3&£nCj a- 
s.t.a>P s.t.<Ta=T 
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X 



n n ) ^ 



vac 



~^ J2 sgn(a,/5)(^[cr] -^[o-^^s]) 



X 



n n 4,.K) 1 $vac. (3.13) 

.xG£\{a} / \fc6j\{j} 



Since this quantity vanishes for all j G JT", we finally find that 

g[cT]=g[a^^f,]. (3.14) 

Since the entire lattice is connected, ()3.14|) ensures that the lowest energy state is unique 

(s) 

in each sector with a fixed S^^i. Therefore $2T is the unique ground state apart from the 
degeneracy for rotational invariance. This completes the proof of Theorem 12.21 

It remains to prove that d states cannot be doubly occupied. Since a states cannot be 
doubly occupied, the ground state can be expanded in the basis states ()3.11|) as 

$Gs= Yl fiE^,E^,JhJl)'^2{E^,E^,J^,J^). (3.15) 

Et^,Ei,J^,Ji^ 
E^nEi=tl) 

Let Uj and Vj be the internal sites which belong to a cell Cj. By using (j3.15|) and the anti- 
commutation relations {c^^.^^, aj,_^} = -(z//2)5^_^xk ^ ^ n Cj], {%,cr,4,T} = ^a,TX['^k e Cj], 
{cvj,cT, 4,r} = -^a,TX[uk ^ Cj] for auy x e£, j,k e J, we get 

C7„TC^„i'^'GS = Y 5Z ^(^T, El, Jt, Ji)sgn(x, y)^2iE^\{x}, Ei\{y}, J^, J^) 

x,yGCj E^,Ei,Jt;,Ji 
s.t.xeEf,y<^Ei 

Ei;nEi=<D 

E E ^(^T, El, Ji, Ji)sgn{x, Uj)^2{EM^}, E^, J^, Ji\{j}) 

x&Cj E^,Ei,.J-^,.Jy 
s.t.xGE^ 
E^-nEi=& 

E E f{Ei,Ei,Ji,Ji)sgn{u„y)^2{Ei,Ei\{y},Ji\{j},Ji) 

yeCj E^,E^,J^,Jy 
s.t.y&Ey 
E,^r\Ei=% 



+ Yl f^^h ^i' "^T, ^i)sgnK, u,)ME^, E^, JMj}, Jl\{j}) (3-16) 

s.t.E-fnEi^=9 

where =f is — for 7 = m and + for 7 = v, and sgn{x,y) depend on the sites x,y and 
Ej, J|, Jj. 

These quantities vanish for all j & J because of (j3.5|) . By adding these two relations 
thus obtained, we find that 

fiE^,Ei, Jt, Ji)sgniu„u,)^2iE^,E^, J^j}, Jl\{j}) 

E-j,E^,J^,Jl 
s.t.EinEi=Q 

+T 5Z E f^^h El, Jt, Ji)sgn(x, y)M/2(i?T\{^}, ^T, Jl) = (3.17) 

x,yeCj E^,E^,J^,J^ 
s.t.x€Ei-,yeEi 
EfnEi=(l) 

The linear independence of implies that each coefficient of ^2 in ()<^.17j) is vanishing. 
This means that 

f{Ei,Ei,J^,Jl)sgn{uj,Uj) = 

E E fiE\,E[,J{,J[)sgnix,y) (3.18) 

x,y&Cj E'^,E[,J'^,J[ 
E'^nE[=iD 
s.t.{E\\{=o},E[\{y},J\,J[} 
={E^,E^,J^\{j},J^\{j}} 

for any Et-,Ei C £ and J|, J| C jT" with E-^ HE^^ = 0. Note that the two electrons occupying 
the same d state in {Ei^,Ei, J-\,Ji} are removed and added to S in the new configuration 
{E\,E[,J\,J[}. 

Now suppose that d state is doubly occupied in r different cells. More precisely we assume 
that there is a configuration {E-^, E^, J|, J^} with /(-E|, E^, J|, 7^ such that | J|n = r. 
We use the relation (j3.18p to this configuration by setting j G J| fl J^. Then one finds that 
the summation in the right-hand side of (j3.18|) is over configurations {E'p E'^, J|, J|} where 
the d state in Cj is not occupied and there are two more occupied a states compared with the 
original configuration. Thus {E'p E'^, J'p J|} has (r — 1) cells with doubly occupied d states 
and at least one cell with the empty d state. Because of ()3.18|) . there is at least one such 
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configuration with non vanishing /(E'j, J|, J|). Repeating the argument r times, one 
concludes that there is at least one configuration {E'p E'^, J|, Jj} with E'^, J|, J|) 7^ 

which has no cells with doubly occupied d states and at least r cells with no d states. But 
this is a contradiction since the maximum possible electron number for such a configuration 
is \S\ + \ J'\ — r < \£\ + \ J\ = Ne. This prove the desired claim that r is always 0. 

4 Discussions 

Let us make two remarks about our model. 

First we discuss the mechanisms that generate ferromagnetism in the present and Tasaki's 
models. In Tasaki's models, electrons in the lowest fiat-band may be regarded (in the basis 
corresponding to the a operators) as almost localized at external sites. Roughly speaking, 
a small overlap of the wave functions at an internal site generates "exchange interaction" 
which leads to ferromagnetism as in Figini(a). In the situation of Theorem 12. the picture 
is almost the same in our models. The electrons are almost localized at external sites and 
overlap at intermediate sites as in FiglHl (b). In the situation of Theorem \2.2\ however, the 
picture is essentially different from that in Tasaki's model. Each electron in the lowest flat 
bands is either almost localized at an external site or localized at a pair of internal sites. 
The basic "exchange interaction" involves three electrons as in the Figini(c). This is why 
the proof of Theorem 12.21 required a new technique. 




Figure 3: Exchange mechanisms in the one dimensional models, (a) An almost localized 
state for electrons in the lowest flat-band (a state), (b) In the case of Theorem 12.11 overlap 
of two states almost localized at neighboring external sites generate ferromagnetism just as 
in Tasaki's flat-band model, (c) In the case of Theorem 12. 2[ the exchange interaction also 
involves a state localized at two internal sites. 

Secondly let us discuss the possibility of further extending Tasaki's construction. A 
natural question is whether one can treat models with three or more internal sites. As for 
results corresponding to Theorem 12. H it is obvious that our proof (and Tasaki's original 
proof) automatically extends to such models. But results corresponding to Theorem 12.21 
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which involves a new exchange mechanism, are much more dehcate. We suspect that a new 
idea is required to cover the general cases. 

Part of the present work was done when the author was at Department of Physics, 
Gakushuin University. I wish to thank Hal Tasaki for introducing me to the problem of 
ferromagnetism and for various useful discussions and comments, Akinori Tanaka for in- 
dispensable comments, suggestions and discussions. I also thank Masahito Ueda for useful 
discussions and continual encouragements. 
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